Collective modes of CP^ Skyrmion crystals in quantum Hall ferromagnets 
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The two-dimensional electron gas (2DEG) in a bilayer quantum Hall system can sustain an in- 
terlayer coherence at filling factor v = 1 even in the absence of tunneling between the layers. This 
system, which can be described as a quantum Hall pseudospin ferromagnet, has low-energy charged 
excitations which may carry textures in real spin or pseudospin. Away from filling factor i/ = 1, a 
finite density of these is present in the ground state of the 2DEG and forms a crystal. Depending on 
the relative size of the various energy scales, such as tunneling (A^as), Zeeman coupling (Az) or 
electrical bias (A;,), these textured crystal states can involve spin, pseudospin, or both intertwined. 
This last case is a "CP'^ Skyrmion Crystal". In this article, we present a comprehensive numerical 
study of the collective excitations of these textured crystals using the Generalized Random-Phase 
Approximation. For the pure spin case, at finite Zeeman coupling the state is a Skyrmion crystal 
with a gapless phonon mode, and a separate Goldstone mode that arises from a broken (7(1) sym- 
metry. At zero Zeeman coupling, we demonstrate that the constituent Skyrmions break up, and the 
resulting state is a meron crystal with four gapless modes. In contrast, a pure pseudospin Skyrme 
crystal at finite tunneling has only the phonon mode. For Asas — > 0, the state evolves into a meron 
crystal and supports an extra gapless (t/(l)) mode in addition to the phonon. For a CP^ Skyrmion 
crystal, we find a U{1) gapless mode in the presence of non-vanishing symmetry-breaking fields 
AsAS, Az, and A;,. In addition, a second mode with a very small gap is present in the spectrum. 
We present dispersion relations for the different low-energy modes of these various crystals as well 
as their physical interpretations. 

PACS numbers: 73.21.-b,73.22.Lp, 73.20.Qt 
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I. INTRODUCTION 

The two-dimensional electron gas (2DEG) in a bilayer 
quantum Hall system has a broken-symmetry ground 
state at filling factor ly — 1 that can be described as 
a pseudospin ferromagnet. This state is characterized by 
finite interlayer coherence even in the absence of tunnel- 
ing. Interlayer coherence is maintained if the interlayer 
separation d is lower than a critical separation which is 
of the order of dc « 1.25^ at zero electrical bias where 
£ = y/ Hc/eB is the magnetic length. This state has been 
extensively studied, both experimentally and theoreti- 
cally (for a review, see Ref. [ll). Its collective excitation^, 
which has been detected experimentally"^, is a pseudospin 
wave mode in which the pseudospins precess around their 
equilibrium orientations. 

In most studies of the interlayer coherent states at 

= 1 or at other filling factors where states such as 
Wigner crystals, Skyrmion crystals or bilayer stripes are 
predicted, it is assumed that the 2DEG is spin polarized 
so that the non-interacting 2DEG can be mapped onto a 
two- level system with bonding and anti-bonding states. 
Recent experiments^!^, however, cast some doubt on the 
validity of this assumption. To explain these recent ex- 
perimental findings, it seems that the ground state should 
be partially spin depolarized. 

A previous study^ (which we refer to as Paper I in the 
rest of this article) presented various zero-temperature 
phases of the 2DEG in a bilayer quantum Hall system 
when the filling factor is varied around v = 1, so that a 



finite density of charged excitations is introduced in the 
ground state. Working in the Hartree-Fock approxima- 
tion, various symmetry breaking fields were present in 
the system: tunnelling (Ag^g), Zeeman (A^) and elec- 
trical bias (Af,), over a range of experimentally relevant 
filling factors (0.6 < < 1.2) and interlayer separations 
(0 < d/£ < 1.0). The phase diagram depends sensitively 
on the interplay of the different symmetry breaking fields 
as well as filling factor and interlayer separation. In some 
regions of the phase diagram, the 2DEG is pseudospin po- 
larized and the charged excitations are Skyrmions with 
real spin texture (i.e., spin-Skyrmions) while in regions 
where the Zeeman coupling dominates, the 2DEG is 
spin-polarized and the charged excitations are Skyrmions 
with pseudospin texture (i.e., pseudospin Skyrmions or 
bimerons.) The phase diagram contains also some re- 
gions where the charged excitations involve intertwined 
spin and pseudospin textures. These complex objects are 
called CP^ Skyrmions''. Spin and pseudospin-Skyrmions 
can be viewed as limiting cases of the more general CP^ 
Skyrmion excitation. 

In the present paper, our main goal is to under- 
stand the collective mode spectrum of a crystal of CP^ 
Skyrmions. To understand the intricate spectrum of ex- 
citations in such a crystal, we find it helpful to begin 
our study by calculations of the collective mode spec- 
tra of pure spin or pseudospin Skyrmion crystals as well 
as that of the homogeneous (or liquid) state at exactly 

For the liquid state, an analysis of the collective ex- 
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citations was presented in Ref. d, but the analysis was 
limited to the behavior at long wavelengths. We provide 
here the full dispersion relations, calculated in the Gen- 
eralized Random-Phase Approximation (GRPA), of the 
three Goldstone modes identified in Ref. 8. In addition 
to pseudospin-wave and spin-wave modes, the CP^ liq- 
uid state also supports a gapless mode involving a spin 
and a pseudospin flip. 

The collective mode spectrum of apure spin-Skyrmion 
crystal was first calculated in Ref'y. In this work we 
present and discuss both two-component cases: the spin 
and pseudospin crystals, and discuss the similarities and 
differences. The spin-Skyrmion has two Goldstone modes 
at finite Zeeman coupling. One is the phonon mode 
which is due to the broken translational symmetry, while 
the second mode is a U {!) phase mode arising from a 
broken orientational symmetry in the spin texture. In- 
deed, the energy of an isolated spin-Skyrmion is invariant 
with respect to a global change in the azimuthal angle 
of the spins. In a crystal, this angle is fixed and this 
symmetry is broken. It is believed^ that this mode is 
responsible for the small NMR relaxation time Ti ob- 
served experimentallyi^ around = 1 at low temper- 
ature. A particularly interesting situation arises if one 
considers the limit of vanishing Zeeman coupling A^: 
each Skyrmion breaks up into two merons (a meron is 
a vortex-like configuration which can be understood as 
half a Skyrmion), and the resulting meron crystal has 
three Goldstone modes in addition to the phonon mode 
which are similar in nature to the U {1) phase mode of 
the spin-Skyrmion crystal. 

For the analogous bilayer situation, where the two lay- 
ers in which the electrons may reside play the role of the 
two allowed spin states, we find that only the phonon 
mode remains gapless in the pseudospin- Skyrmion crys- 
tal when the tunneling and layer separation are both 
non- vanishing. The absence of a gapless mode analogous 
to the t/(l) mode of the spin-Skyrmion crystal occurs 
because the inter-layer interaction is different than the 
intra-layer interaction, breaking the U{1) symmetry. For 
AsAS — ^ 0, the state deforms into a pseudospin meron 
crystal, which has two (rather than four) gapless modes 
because of the explicit symmetry-breaking in the inter- 
action. These modes may be identified as a phonon, and 
a Goldstone mode due to the spontaneous interlayer co- 
herence in the limit of zero tunneling and the broken 
orientational symmetry in the in-plane pseudospin tex- 
ture. 

For the CP^ Skyrmion crystal, we show that, when all 
symmetry breaking fields Asas, A^, and Af, are finite, 
the crystal has, in addition to the phonon mode, one 
gapless phase U {!) mode arising from the invariancc of 
the Hamiltonian with respect to a rotation of all the spins 
in the two layers by the same angle. We find also that 
since the CP^ Skyrmion crystal occurs only at very small 
value of the tunneling parameter, a second phase mode 
with a very small gap is present in the spectrum. Our 
CP^ Skyrmion crystal is characterized by a set of eight 



low-energy modes that offer a distinct signature of this 
entangled spin and pseudospin state, much different from 
that of a pure spin or pseudospin Skyrmion crystal. 

Our paper is organized in the following way. We intro- 
duce our description of the crystal phases in the bilayer 
quantum Hall system in Sec. II and the formalism we 
use to compute the dispersion relations of the collective 
modes in the GRPA in Sec. III. We study the collec- 
tive mode spectrum of the liquid state for v = 1 'm Sec. 
IV, that of the pure spin-Skyrmion crystal in Sec. V, 
and that of the pure pseudospin crystal in Sec. VI. We 
are then ready in Sec. VII to consider the more general 
case of a CP^ Skyrmion crystal. We study the excita- 
tion spectrum of this crystal as a function of interlayer 
separation, tunnel coupling and electrical bias. Finally, 
we conclude in Sec. VIII. 



II. DESCRIPTION OF THE CRYSTAL STATES 

We consider a 2DEG in a double-quantum- well system 
(DQWS) in a quantizing magnetic field B = Hz, taking 
into account the possibility that an electrical bias is ap- 
plied to the system. We define this bias as A^ = En — El 
where E^i^l) are the electric subband energies in each well 
(right and left) in the absence of magnetic field and tun- 
neling. For simplicity, we make a narrow well approxima- 
tion, i.e. we assume that the width h of the wells is small 
{h « ^) and treat interlayer hopping in a tight-binding 
approximation. The single-particle problem is then char- 
acterized by the total filling factor z/, the separation d 
(from center to center) between the wells, the tunneling 
gap As^s, the electrical bias A{, and the Zeeman energy 
A^ ~ g* ilbB where g* is the effective gyromagnetic fac- 
tor and /is is the Bohr magneton. 

We describe the phases of the electrons in the DQWS 
by the set of average fields |^Pi^/(q)^| where p"^'^(q) 
is an operator that we define^ as 

Pi,j \H) — ^ 2_^e C^^^^x'^j,l3,X-qyP- 

(1) 

In this equation, — S/2'itI'^ is the Landau level de- 
generacy with S the area of the 2DEG. The indices 
i,j — R,L and a,/? = ±1 are respectively well and 
spin indices. We use the Landau gauge with vector po- 
tential A = (0,i3x,0). The operator c\ ^ ^ creates an 
electron in state {i,a,X) in the lowest Landau level. 
We work in the strong magnetic field limit where the 
Hilbert space is restricted to the lowest Landau level 
only and make the usual approximation of neglecting 
Landau level mixing. For a crystal state, only the fields 

Pifi^ ~ 7^ 0' where G is a reciprocal lattice vec- 
tor of the crystal structure. The calculation of the order 
parameters (^p"f{(\ — G)^ in the Hartree-Fock approxi- 
mation for a crystal of CP^ Skyrmions was explained in 
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Paper I^. 

The Hartree-Fock Hamiltonian of a general crystal 
state is given by 



Hhf = iV0 5]s„,.p"f (0) 



(2) 



A 



SAS 



E E E (G) pf,f (g) 

Q,/3 i,j G 

-^^ E E E (G) 4"(G), 

Q,/3 i,j G 



with the renormalized single-particle energies given by 



El,, 



2 

2 
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Vi) V2 
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e 



(3) 
(4) 



and the Hartree and Fock intrawcU and interwell inter- 
actions defined by 



H^,^ [q) 


= H{q) 




= H{q) 


(?) 


= X{q) 




= X{q) 











dye-^ /Vo (g^y) ,(7) 
(8) 

+c>o 



Note that to take into account a uniform neutralizing 
positive background, one sets H {{)) = and H (0) = 

For the full four-level system, an electronic state may 
be specified by the four-component spinor 



cx 




(9) 



At this point, it is convenient to redefine the opera- 



tors p^f (q) introduced above as Pa,6(q) with the indices 



a, 5 taking the four values a,b = 1,2,3,4 and referring 
to the states (i?, , (i?, — ), (L, (L, — ) in this order. 
The layer index can be mapped into a pseudospin index 
ibpwith the correspondance R —^ +p and L ^ —p. The 
density, spin density, and pseudospin density operators 
p (q) , S (q) and P (q) and the nine operators R^^y (q) 
(with u,v = X, y, z) are relatedi^ to the four-component 



spinor of Eq. ([9]) by the equationsi^ 

^(q) = ^E^""^"'' 



/2 c\cx-,^p, (10) 



X 



P, 



-iq^X+iq^qyfi /2 t spin 



^X^u CX-q^P, 

(11) 



(n\ — p-iq^X+iq^qyi'^ /2 t ppin 



X 



(12) 



P (ci) — ^ 'Sr^ -iq^X+iq^qyt' /2 J spin ppin 
nu^vy'A) ~ rfjsj 2-^^ '^X^u ^X-qyl^: 



X 



(13) 

where the 4x4 matrices r^^*" and r^^**" are defined by 

cr„ 



spin 



CT„ y ' 

(with cr„ a Pauli matrix) and by 



(14) 



^pp^n = ( { y ^PP^n ^ f 0^ 1 ^ (^5) 



/ / ' 



and 



/ 

-/ 



(16) 



where / is the 2x2 unit matrix. 

One may show^ that the order parameters obey the 
four sum rules 



(17) 



G b 



III. CALCULATION OF THE RESPONSE 
FUNCTIONS 

We define the two-particle Matsubara Green's func- 
tions as 

Xa,b.c,d (q, q'; t) = -N^ {TrSpa,b (q,T) Spc.d (-q', 0)) , 

(18) 

where Tr is the imaginary time ordering operator. In a 
crystal, these functions are non zero only for wavectors 
q = k -t- G and q'= k -I- G' where k is a vector in the first 
Brillouin zone of the lattice. The collective mode spec- 
trum is found by tracking the poles of the response func- 
tions Xa.b,c,d (k -|- G, k -|- G; i^n —^i-u + iS) when k varies 
in the Brillouin zone. 

Let us denote by a; , the layer and spin indices in a 
and similarly for the other indices 6, c, d. The equations 
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of motion for the two-particle Matsubara Green's func- 
tions in the Hartree-Fock approximation (i.e. Xa b c d) 
given by 



(19) 



= -N^fi {[PaM (q,r) , pc^d (-q', 0)]) 6 (r) 

-N^ {Tr [Hhf - M^, Pa^b (q,T)] , 5p,,d (-q', 0)) , 

with Hh f given in Eq. ^ . Evaluating the commutators 
and Fourier transforming with respect to the imaginary 
time T, we get the lengthy equation 



n + iEa -Eb)] x°,6,c,d(q:q';f^r, 

{Pa4 (q - q')) f^bxTq.q' 
-fl{Pc,b (q- q'))'5a.d7q,q' 



(20) 



^SAsSbi,RXa,(b,,L),c,d i'ii ^ ^ri 



-^AsAsSb,,LXa,{bs,R),c,d (q> q'; ^«) 

-EE <a(q-q")7q,q"x%,c,<i(q",q';f^«) 

a' q'Vq 

+ E E t^i^,aq"q")7WA,c,.(q",q';f^n) 

a' q'Vq 

+ E E (q - q") 7q,q"x?' Acrf (q"' q'; ^n) 

b' q" 

-EE ^^'^b (q - q") 7*,q"X°,a'c,d (q", q'; ^n) , 

a' q" 

where 7q.q' — e*^''^'')'^^ and we have defined the 
mean-field Hartree and Fock potentials 

C/5(q) = (q) (p"-" (q)) , (21) 

C^M(q) = X<,,^,,(q)(p-'^(q)), (22) 

with the Hartree and Fock interactions defined in Eqs. 

dMHl). 

We get the response functions in the GRPA by sum- 
ming the ladder and bubble diagramsii. The result is 
the integral equation 

XQ,ft,c,<i(q,q';^^n) (23) 
= Xa,b,c,d (q, q';f^n) 

+ ^ E ^-Ae.e (q, q"; ^n) H^.g (q") Xg,g.cM (q", q'; ^n) 
q" 

I E ^Ib.ej (q, q"; ^ej (q") xf,e,c,d (q", q'; ao . 



n 



q" 



We now write Eq. ([23]) in a matrix form more 
appropriate for numerical analysis. We first rede- 
fine Xa.b.c,d as the matrix Xi.j with the indices i,j 



taking the values from 1 to 16 with the order: 
(1, 1), (1,2), (1,3), (1,4), (2, 1), etc. for i, the row index, 
and the order: (1, 1), (2, 1), (3, 1), (4, 1), (1, 2), etc. for 
j, the column index. It is to be noted that all ele- 
ments Xi,j of tl^c matrix x well as of the other ma- 
trices I,A,B,H,X that we define below are themselves 
matrices in the reciprocal lattice vectors G,G'. In the 
end, we have to diagonalize a matrix with dimensions 
ISn^jLy X WriRLv where urlv is the number of recip- 
rocal lattice vectors G that we keep in the calculation. 
We typically work with riRi^y « 317 (in the calculation 
of the order parameters, however, we keep more than 600 
reciprocal lattice vectors). We then write the index i as 
i = 4 (i — l) + i where i,i — 1^ 2, 3, 4 and similarly for j . 
The unit matrix is written as 



I = %'%J'^q,q'- 
We define the following matrices: 

s»j(q,q') (pij(q-q'))7q.q'% 
-(pj,T(q-q'))7q,q'%, 



(24) 



(25) 



(26) 



F^,, (q, q') = - [H+ (q, q') - i?f (q, q') j <5., (27) 
+^rt.(q,q')5jj-^(q,q')%, 



^,3 



rr I l\ ^SAS ^ r 

7^J (q,q ) = — :^^7-d 



2 ■'ij"|j-j|,2"q:q ' 

-^%^|I-7|.2'5q.q' 



(28) 



with the definitions (with fc, Z = 1,2,3,4) 

i?^ (q,q') =E^i^.'c(q-q')e^<•'^•'")■^^^/^ (29) 



and 



(q-q") = t/,^,(q-q')e±<''^'>")-^^^/^. 



We also define the interactions 
V.^, (q, q') - \h^^ (q) - X-r- (q) 



(30) 



With all these definitions, we can write the equation of 
motion for the response functions in the matrix form: 



[iMlnl - E -T - F - BV]x^fiB. (31) 

To get the weight of each pole, we diagonalize the matrix 
M = E + T + F + BV hy writing MU = UD where U 
is the matrix containing the eigenvectors and D is the 
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diagonal matrix of the eigenvalues Hujj (k). The response 
function x can be written as 



OJ - 



i6 



(k) 



(32) 



where the indices n, m, j, k — 1,2, Wuhlv Using Eqs. 
P0m3p . it is easy to get from x the response functions for 
the density, spin, pseudospin and the R fields. Because 
we have access to the matrix of eigenvectors U, we can 
also animate the motion of the density, spin, pseudospin, 
etc. in a given mode. This will be very helpful for the 
interpretation of the different modes. 



IV. UNIFORM COHERENT STATE AT i/ = 1 

In this section, we apply our formalism to the study 
of the spin-polarized uniform coherent state (UCS) at 
1/ — 1. For effectively spinless electrons, the behavior of 
the collective mode (i.e. pseudospin wave) with bias and 
interlayer separation was studied before in Ref. [3. One 
may also consider— the limit of an S'[/(4) symmetric 
Hamiltonian, in which the interlayer separation d, Az, 
and As^is are all set to zero. In this case three Gold- 
stone modes were predicted in the absence of symmetry- 
breaking fields^ using a gradient expansion, but quanti- 
tative dispersion relations for these modes that takes into 
account the spin, pseudospin and density fields are dif- 
ficult to obtain in this approach. Our GRPA formalism 
produces these three modes, provides analytical results 
for the dispersions in some limiting cases, and allows nu- 
merical results to be obtained in the general case from 
the microscopic Hamiltonian. 

In UCS, the only non-zero order parameters are 



{pi.i m 

(P3,3 (0)) 
(Pl,3 (0)) 



i(P3,i(o))r 



(33) 



a, 



(we can choose a real without lost of generality). For a 
balanced DQWS, i/i — while in the presence of a bias, 
vi , V3 and a must be found by minimizing the Hartree- 
Fock energy. We have discussed the effect of a bias on 
the UCS in Paper I^. In the absence of tunneling, we 
obtain a simple relation between charge imbalance and 
bias, 



-A. 



a = 1^3 - i/i 



x(o)-x(o)-f(i) 



(34) 



In the UCS, the matrices in Eq. ([3T|) have dimen- 
sions 16 X 16 and the system of equations split into three 
uncoupled systems of equations corresponding to the 16 
excitations represented schematically in Fig. 1 (transi- 
tions 1 to 4 are intra- level excitations with zero energy) . 
The excitations represented with the dashed lines in this 
figure are just the reverse of those represented with the 
full lines. In this figure, we have illustrated the situation 



where Ah = and As^g < Az so that the first two lev- 
els are the symmetric (S) and antisymmetric (AS) states 
of the DQWS with real spin up. Because only the S, + 
level is filled in the ground state, our calculation gives 
three dispersive modes corresponding to the transitions 
from the filled level S, + numbered 5, 11 and 13 in Fig. 
1. We will not discuss the other non dispersive modes 
as they have zero weight in the response functions at 
zero temperature. Note that, at finite bias, the bond- 
ing and antibonding states (with single-particle energies 
A'g^g + A^) replace the S and AS states in the dia- 
grams of Fig. 1. 
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FIG. 1: The 16 possible transitions between the single-particle 
energy levels in the UCS is divided in the GRPA into 3 un- 
coupled subsystems involving 8,4 and 4 transitions. 

Transition 5 corresponds to a pseudospin-wave mode 
first found in Ref. 2 in which the pseudospins precess 
around their groundstate orientation. Solving Eq. (|3ip . 
we find the general dispersion relation 



h^Loliq) = 4a (g) 5(g) 



A, + a(X(0)-X(g)-^(^ 



(35) 
-I 2 



where 



aiq) 



^5A5 

2 

AsAS 



X{Q)-X{q) 



(36) 



-fa 



X{Q)-X{q) + H{q)-H{q) . (37) 



Using Eq. (p4|) . it is easy to show that this mode is gap- 
less if As^s = even with finite bias. This is true as 
long as the interlayer coherence is non zero i.e. a ^ 0. 
At q = 0, the pseudospins execute small oscillations in 
X ^ y plane if the bias is zero while, at finite bias and 
q — Q, they execute small oscillations in a plane where 
{Pz) — est. In both cases, the energy cost is zero because 
the Hamiltonian of Eq. ([2]) is invariant with respect to a 
rotation in a plane of constant (Pz) if As as = 0. At fi- 
nite tunneling, the pseudospin-wave mode is gapped with 
^SAS' frequency at g = 0, larger than Asas- 

The dispersion relations of the collective modes num- 
bered 9,11,13,15 are found by solving detB — with 
the matrix B given by 
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B 



lyiX (q) 
2 

- aX (q) 




As 



/ 



Cl 



- >^lX (q) - 
0^ 

f-aX (g) 



A, 



C3 





v^X {q) + Ab 



/ 



\ 



AsAS 

2 



(38) 



C3 - J/gX (g) / 



In this equation, we have defined the constants 

Cl = Az + viX (0) , (39) 



C3 = Az + J^3^ (0) : 



/ = ^+«^(0), 



Ab = A6 + (1/3-1^1) 



(40) 



(41) 



(42) 



For zero bias (in which case a = 1/2), we can solve ana- 
lytically to find the two dispersive modes 



fiwii {q) ^ Az + /S.SAS 

1 

+ 2 

fiwis [q] = Az 
1 

+ 2 



X (0) + X (0) - X (g) - X (g) 



X (0) + X (0) - X (g) - X (g) 



(43) 



(44) 



The frequency 0^13 (g) corresponds to the spin-wave (or 
Zeeman) mode in which the spins precess around their 
equilibrium orientations while mode 11 involves a spin 
and pseudospin flip. At zero bias and Ag^is = 0, these 
two modes are degenerate at g = i.e. fkjon (0) = 
?iwi3 (0) = Az. 

Fig. 2 shows the evolution of the dispersion of the 
three modes 5,11,13 with bias at d/i = 1,A5A5 = 
0, Az/(eVK£) = 0.2 for (a) A^/ (cVk^) = 0, (b) 
Ab/ {e^/Kt) = 0.1, and (c) Afc/ {e^/Kt) = 0.4. For case 
(d), Ab/ {e'^/ni) = 0.4 and Ag^s/ {e^ / ntj = 0.02. Our 
numerical method also produces the frequency of the 
non-dispersive transitions 7,9, 15. We are not interested 
in these transitions as they have zero weigth in the re- 
sponse functions at zero temprature. 

From Fig. 2, we see that there is one pseudospin mode 
and two spin-wave modes. We find that all three disper- 
sions are affected by the bias. Figure 2(a) shows a roton 
minimum A;,/ (c^/k^) — 0.4 increasing with interlayer 
separation indicating the instability of the UCS. Figures 
(b) and (c), show that an increase in the bias removes 
this roton mimimum. 

The value A;,/ {e^ /k,£) = 0.2 in Fig. 2(b) corresponds 
to the filling factors vi = 0.38, 1^3 = 0.62 while the value 
Afc/ [e^/Ktj = 0.4023 in Fig. 2(c) corresponds to i^i = 



0.00008,1/3 = 0.99992. Thus, the bias in Fig. 2(c) is 
just below the critical value above which the charge is 
completely transferred to the left well. It is easy to get 
analytical results for the asymptotic values (g — > 00) of 
the modes in this case (i.e. the transport gaps). We 
need only solve Eq. ([55]) with t = Q,iyji = and take 
X (q) ,X (q) — > 0. We find a spin-wave (SW) and a 
spin-wave with pseudospin flid (SWPF) with dispersions 
fioJswPF (g 00) = Ab + Az - j{e^/Ke) + X (0) and 
hujsw (9 ~^ 00) = Az+X (0) . In this limit, the S and AS 
symmetric states have evolved into the L and R states. 
The mode SW corresponds to a spin flip transition in the 
same well and its energy is fkosw il 00) = A^ -1--'^' (0) 
since X (0) is the energy required to remove an electron 
from the filled L, | state. This mode is the spin-wave 
mode in a single quantum well system. Its dispersion 
relation at finite g is given by 



Hiusw (q) = Az+X{0)-X (q) , 



(45) 



The mode SWPF involves a spin and a pseudospin flip 
and so a transition of an electron from the L to the R 
well. When the charge is completely transferred into 
the left well and we take into account the positive neu- 
tralizing backgrounds, there is an electric field oriented 
from the R well to the L well. There is thus a gain in 
energy j{e'^ / nd.) in transferring the charge from the L 
to the R well. The total energy to remove an electron 
from the filled L, t state and create an infinitely sepa- 
rated electron-hole pair with spin and pseudospin flip is 
thus hujswPF (q ^ 00) = Ab + Az - j{e^ I nt) + X (0). 
With the parameters of Fig. 2, we find fkosw (q ~* 00) = 
1.25{e^/Ki) and fujswPF (g -> 00) = OM{e^/Ki). The 
values are in good agreement with the numerical result 
plotted in Fig. 2(c). In Fig. 2(d), we have added a small 
tunnel couphng AsAsHe^/n^) — 0.02 that fifts the de- 
generacy of the two modes gapped at the Zeeman energy. 
At g = 0, the ordering of these two modes is just the op- 
posite than at g — > 00. The Zeeman mode is gapped at 
fkosw (0) = Az while the pseudospin-|-spin flip mode is 
gapped at ftioswPF (0) ^ Az + Asas- 



V. SPIN-SKYRMION CRYSTALS 

If the electrical bias is sufficiently strong for all the 
charges to go into the left well, the system becomes effec- 
tively a single quantum well system (SQWS). The index i 
then takes the values i = 3, 4 only. It is well known that. 
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FIG. 2: Dispersion relations of the collective excitations in 
the UCS at = 1 as computed in the GRPA. Parame- 
ters are dji = 1, Az/ {e^ j nt) = 0.2 and (a) Ab/ (f? j ni) = 
0,AsAS = 0, (b) A6/(eV«:^) = 0.1,AsAS = 0, (c) 
Ab/ (eV^^) = 0.4023, Asas = 0, and (d) At/ {(? j nt) = 0.4, 
AsAs/ (e^/K^) = 0.02. The straight lines give the frequency 
of the non-dispersive modes. The modes 5, 11, 13 correspond 
to the dispersions cus (q) (q) ,0^13 (q) derived in the text. 
The abbreviations SW and SWPF refer to the spin-wave and 
spin-wave with pseudospin-flip modes (see text). 



in this case, the ground state at = 1 is a spin ferromag- 
net. At small Az, this state supports topological spin- 
texture excitations known as SkyrmionsiSiii^. Away from 
I' = 1, a finite density Ug = \u — 1\ /2Tr£'^ of Skyrmions 
is included in the the ground state and is expected to 
form a Skyrme crystal— . The Hartree-Fock ground state 
energy of the Skyrme crystal in a SQWS is, using Eq. 



E 



HF 



N 



(46) 



+^E^2(G)i(p(G))r 

G 

-iE^(G)i(s(G))r 



G 



where we have defined 



T2 (G) = 2F(G)-X(G) 



(47) 



A Skyrme crystal has a noncoUinear magnetic order. 
A single Skyrmion spin texture of such a crystal has its 
spins aligned with the Zeeman field at infinity, reversed 
at the center of the Skyrmion, and has nonzero XY 
spin components at intermediate distances which have 
a vortex-like configuration (see Fig. 3(a)). The classical 
(or quantum mean field) energy of a Skyrmion is indepen- 
dent of the angle (p which defines the global orientation 
of the XY spin component. This independence gives an 
extra U {!) degree of freedom for a single Skyrmion. In 
a crystal of Skyrmions, this U (1) symmetry is sponta- 
neously broken and the crystal has an extra Goldstone 
mode. We refer to this mode as the XY mode because it 
corresponds to an oscillation of the perpendicular com- 
ponent of the spins in the x — y plane with (Sz) fixed. 
From Eq. (|46p . we see that such a motion costs no energy 
provided the density {p{G)) is not changed. That this 
is the case can be seen from the fact that the relation 
between the topological charge and the spin texture^- is 
given at 1/ = 1 by (the sum over indices is implicit) 

S {n (r)) = -^SabcSa (r) £ijdiSb (r) djSc (r) , (48) 
an 

where and Sabc are antisymmetric tensors, with i,j = 
X, y and a,b,c = x, y, z, and s (r) is a classical field with 
unit modulus representing the spins. If we write a general 
spin texture as 



Sx (r) = sin [0 (r)] cos [ip (r)] , 
Sy (r) = sin [9 (r)] sin [(p (r)] , 
Sz (r) = cos [9 (r)] , 

then the induced density takes the simple form 
1 



6n (r) 



An 



sin [61 (r)] [V9 (r) x V(/j (r)] -2. 



(49) 
(50) 
(51) 



(52) 
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Clearly, the induced charge density in the XY mode is 
zero because Vip (r) is unchanged if we rotate all the spins 
by the same angle (i.e. (p{r) ^ (p (r) + est). 

The phase diagram of the 2DEG around 1/ — 1 has 
been studied in the HFA in Refs. [g UlTlfigl . In a large 
portion of the i' — Az phase space, the Hartree-Fock 
ground state is a square lattice with two Skyrmions of 
opposite global phase tp per unit cell as shown in Fig. 
3(a). This configuration is called SLA (Square Lattice 
Antiferromagnet) to account for the phases (p = and 

= TT of the two Skyrmions in the unit cell. As the Zee- 
man energy is increased, our Hartree-Fock approximation 
(HFA)^- shows that there is a transition to a triangular 
lattice of Skyrmions with a phase difference of (/? = 27r/3 
between the three Skyrmions in a unit cell and then to a 
triangular Wigner crystal of spin down electrons with no 
spin texture. 

Note that the size of the individual Skyrmions in the 
Skyrme crystal is determined in part by the Zeeman 
energy Az- For large Az, the tipping of spins into 
the minority direction that occurs in the spin texture 
of a skyrmion becomes energetically unfavorable, and 
Skyrmions become relatively small. For decreasing Az, 
the Skyrmions become larger until their size is limited 
by inter-Skyrmion interactions. For small enough Zee- 
man coupling the system undergoes a transition into a 
meron crystal^^ with four merons per unit cell as shown 
in Fig. 3(b). 

A meron is most easily visualized by rotating the spin 
axes around the y axis so that the z axis becomes the x 
axis and the x axis becomes the —z axis. In this basis a 
meron has its spin down or up at the origin and lies in the 
XY plane at infinity with a vortex-like structure. The 
vorticity can be positive or negative so that, combined 
with the up or down spins at the origin, there are four 
possible types of merons. In the rotated basis, one can 
recognize the Skyrmion as a bound meron-antimeron pair 
(i.e., a bimeron). Hartree-Fock calculations^Si shows that 
the meron crystal is stable only at filling factors ly > 1.38. 
For smaller filling factors, we find, as suggested in Ref. 
[21I, that Skyrmions always appear in pairs to form a tri- 
angular Wigner lattice. There is an attractive interaction 
between Skyrmions with opposite phases which prevails 
over the Coulomb repulsion at large separations between 
the Skyrmions. Thus, as the Zeeman energy is reduced 
and for < v < 1.38 where Vc ~ 1.06, there is a transi- 
tion from the Skyrme crystal to a crystal of biskyrmions. 
In Ref. [2ll . the minimal energy is found for zero sepa- 
ration between the pairs of Skyrmions in a biskyrmion 
leading to a Skyrmion with topological charge Q = 2. 
Within the HFA, there is a finite separation between the 
two Skyrmions in each pair. 

The low-energy collective modes of the Skyrme^ and 
meron crystals computed in the GRPA are shown in Fig. 
4. The dispersions are shown for wavevector k along the 
path T — M — X — Y around the irreducible Brillouin zone 
as shown in the inset of Fig. 4(b). The low-energy part 
of the spectrum consists of two gapless (Goldstone) and 
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FIG. 3: (Color online) Spin texture in the x — y plane (vec- 
tors) and density in a unit cell of a Skyrme lattice in (a) 
the SLA and (b) the meron configurations at z/ = 1.30 with 
Az/ {e^/K,l) = 0.01 and Az/ {e^/Kt) = respectively. In (a), 
Sz = — I at the position of the Skyrmions and each Skyrmion 
carries charge q — — e. In (b), Sz — a.t the positions of the 
four merons at the corner of the unit cell and at the position 
of the meron at the center while Sz ~ +\ at the positions 
of the other merons. Each merons carries a charge g = — |. 
Contours with dark (red) regions indicate a high density. 



two gapped modes. One of the Goldstone mode is the 
phonon mode which is due to the spontaneously broken 
translational symmetry in the crystal. Its dispersion is 
Lo ~ k^/"^ which is similar to that of a Wigner crystal in 
a magnetic field. The second Goldstone mode is the XY 
mode introduced above. This gapless mode is believed 
to be responsibleSiSS, for the small spin-lattice relaxation 
time Ti measured in NMR experiments^" on both sides 
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oiv=l. 

Of the two gapped branches, one has a gap that varies 
with Az and v and the other is gapped at exactly the 
Zeeman energy for aU v as required by Larmor's theorem. 
This last mode is obviously a spin wave mode where the 
spins precess around their local orientation. We refer to 
this mode as the Zeeman mode. In the first mode, the 
motion of the two Skyrmions with opposite global phases 
ip in a, unit cell (as seen in our animations^) are out of 
phase. We identify this mode as an optical phonon mode. 

At V = 1.2, the two gapped branches are degenerate 
at k = 0. For v > 1.2, the spin- wave branch is connected 
to the phonon branch at point M as in Fig. 4(a) while 
for V < 1.2, it is connected to the XY spin mode. The 
gap in both branches goes to zero as ^ so that, in 
the meron crystal, all four branches are gapless as seen 
in Fig. 4(b). 

We can understand the origin of the three gapless 
modes (in addition to the phonon mode) in a spin-meron 
crystal when Az = by the fact that the energy of such a 
spin texture is invariant with respect to a rotation of the 
spins around three orthogonal directions. Because these 
three symmetries are spontaneously broken in the meron 
crystal, we get three Goldstone modes corresponding, in 
our case, to oscillations in the x — y,y ~ z and x — z 
planes. This interpretation is confirmed by our anima- 
tions of these modes^^. There is no motion of the charge 
in any of these three modes, the optical phonon mode 
present when Az ^ 0, i.e. in the Skyrme crystal, is pre- 
sumably to be found amongst the higher-energy modes 
not considered in our study. As — > 0, the out of 
phase motion of the Skyrmions in the optical phonon 
mode gradually stops until it disappears completely in 
the meron crystal. 

In Fig. 4 and other similar figures, we plot the poles 
with the two (six for the CP^ case below) biggest weights 
in the response functions indicated in the legend. The 
function Xxx stands for xs^.S^, Xyy stands for xSy,Sy etc. 
For very small wavevector k, this gives an indication of 
the nature of the mode. The phonon mode, for example, 
is strongest in the density response Xn,n while the XY 
mode is strongest in Xzz- 



VI. PSEUDOSPIN-SKYRMION CRYSTALS 

If the spin degree of freedom is neglected, then the 
2DEG in a DQWS can be described by its pseudospin 
and charge orders only. In the absence of bias, the ground 
state at = 1 is a pseudospin ferromagnet^^. Yoy nonzero 
interlayer separation d, all pseudospins are forced to lie 
in the XY plane in order to minimise the capacitive en- 
ergy (the term Jz,2 (0) (0)) |^ in Eq. ^ below). 
The charged excitations of this system are pseudospin- 
Skyrmions or bimerons^. A finite density of bimerons is 
included in the ground state when v ^ \. Once again, 
these topological excitations should form a crystal at zero 
temperature. The Hartree-Fock ground-state energy is 
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FIG. 4: (Color online) Dispersion relation of the collective 
modes of the Skyrme crystal for a (a) SLA configuration with 
v — 1.30 and Az/ [e^ j kI) — 0.008 and (b) a meron config- 
uration similar to that of Fig. 3(b) but with v = 1.38 and 
Az/ {e^/K.e) = 0. The dispersion is computed along the path 
F — Af — X — F in the Brillouin zone as shown by the inset in 
Fig. 4(b). The symbols indicate from what response function 
each frequency comes. The softening of the modes seen in (b) 
is due to the meron crystal becoming unstable aX v = 1.36 
with respect to a biskyrmion crystal. 
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written, in the pseudospin language, as 
Ehf '^SAS 



N 



(53) 



where 



and 



+ iET^2 (G)|(p(G))|^ 

G 
G 
G 

J,,2 (G) = if (G) - i? (G) - X (G) ^ 
Pi (G) = P,(G)x + Py(G)y. 



(54) 



(55) 



The interlayer coherence (i.e. (P_l (q)) ^ 0) is lost at 
some critical interlayer separation which is of the order 
of d/£ « 1.2 at zero bias and tunneling. In the coherent 
regime, the phase diagram in the d/£ — Agj^g — v space 
is very rich and has been only partially studied^^^. In 
analogy with the phase diagram of the Skyrme crystal 
in a SQWS, we find, in some regions of the parameter 
space where t^sAS 7^ 0, a solution consisting of a square 
lattice with two bimerons of opposite global phases and 
meron-antimcron orientation reversed per unit cell. We 
call this phase an SLA bimeron crystal. An example of 
such crystal is given in Fig. 5. Near zero tunneling, we 
get a crystal of merons with four merons per unit cell in 
exactly the same SLA configuration as that of the spin- 
meron crystal shown in Fig. 3(b). 

The collective mode spectrum of a bimeron crystal is 
shown in Fig. 6 where again we have plotted the four low- 
est branches only. There is one gapless and three gapped 
modes. The gapless mode is the phonon mode, which 
appears as a pole of and has a a; ^ k^^"^ dispersion 
as in the Skyrme crystal. The next lowest mode has a 
very small gap that closes as the interlayer separation 
d ^ even when tunneling is finite. In this gapless limit 
it has a linear dispersion. This behavior is easily under- 
stood: when d = and Asas = 0, a bimeron crystal 
can be mapped onto a Skyrme crystal with zero Zeeman 
coupling. (A bimeron is obtained by rotating the spin 
texture of a Skyrmion by 90 degrees around the spin x 
or y axis.) If one rotates around the y axis, the result- 
ing bimeron has pseudospins lying along the x direction 
at infinity, with the two merons aligned along the x axis. 
The tunneling term AsasSx for the bimeron crystal then 
plays the role of the Zeeman term AzSz in the Skyrme 
crystal. It follows that, for d — 0, the collective modes 
of the bimeron crystal are identical to that of Skyrme 
crystal with Az = Asas- However, the specific modes 
appear in different response functions. When d = but 
tunneling is finite, Eq. (j53p becomes 
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FIG. 5: (Color online) Pseudospin texture in the x — y 
plane (vectors) and value of Sz (contours) in a unit cell of 
a bimeron crystal in a SLA configuration with v = 1.20 and 
Asas I {(? /i-^l) = 0.02. The full black circles indicate the po- 
sition of the maximum in the bimeron density. Dark (red) 
contours correspond to Sz ~ +\ while gray (blue) contours 
are for Sz ~ — ^. Note that the pseudospin texture for the 
bimeron at the center has a global phase rotation of with re- 
spect to the pseudospin texture of the bimerons at the four 
corners of the unit cell. 



Ehf 
N 



^SAS 



(56) 



G 

i5]x(G)|(P(G)}^ 



G 



From Eq. (j56p . we see that the pseudospins can rotate 
freely in the y — z plane but not in the x — y plane if 
^SAS 7^ 0. The XY Goldstone mode of the Skyrme crys- 
tal becomes a Y Z Goldstone mode in the bimeron crystal 
at d = 0. At finite interlayer separation H (G) ^ H (G) 
and X (G) ^ X (G) and a rotation in the y — z plane 
can no longer be performed without energy cost. The 
Y Z mode is thus gapped at finite interlayer separation. 

We can explain the smallnest of this gap in the fol- 
lowing way. A global rotation of the pseudospins in the 
y — z plane leads to a solid rotation of the bimerons (i.e. 
a rotation of the pairs of merons around the filled black 
circles in Fig. 5). At d = 0, the bimeron density is cir- 
cularly symmetrio2^ and so does not change under a ro- 
tation of the bimeron. At finite d, however, the bimeron 
density is no longer circularly symmetric. A rotation of 
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an isolated bimeron would not change its energy but, in a 
crystal, these charged objects interact with one another 
via the Coulomb interaction, so there is an energy cost 
for such rotations. The energy cost, however, is small 
because it comes from high order multipole interactions. 
(The rotation of the bimerons can be clearly seen in the 
animations'^.) 

At d / but A — (in which case we get a meron 
crystal) , a rotation of the pseudospins in the x — y plane 
without an energy cost is possible. This in this limit a 
meron crystal has a gapless XY mode (see Fig. 7) just 
as a Skyrme crystal. 
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FIG. 6: (Color online) Collective mode spectrum of a bimeron 
crystal in the SLA configuration of Fig. 5 with v = 
1.2, AsAs/ {e^/K,l) = 0.02 and d/£ = 0.1. The path in the 
Brillouin zone is as indicated in the inset of Fig. 4(b) and the 
symbols are as indicated in the legend of that same figure. 

The interlayer-coherent homogeneous phase v = 1 
can sustairtS. a pseudospin-wave mode in which the pseu- 
dospins execute small precession around their local equi- 
librium direction (the x axis for Ag^g ^ 0). The 
pseudospin-wave mode is equivalent to the spin-wave (or 
Zeeman) mode in the Skyrme crystal. We discussed 
this mode before in Sec. IV. Its dispersion is given by 
Eq. ([55]) . For /S.sas = 0, the dispersion (g) is lin- 
ear in q for small wavevector. This mode has been de- 
tected experimentally'^. For As^is 7^ and d = 0, the 
pseudospin-wave mode is gapped at ojpsw (0) = As^is 
as can be seen from Eq. (|35p . At finite d, that gap is 
increased. In Fig. 6, d is small and we see two modes 
gapped at a value near Ag^g. If we replace {pi^j, (0)) , vr 
and vl by their values in the crystal phases in Eq. ((35)) . 
we get a very good estimate of the pseudospin-wave gap 
uj^ (0) in the bimeron crystal. In the example of Fig. 
6, d/£ = 0.1, (pi,3 (0)) = 0.106 and = = 0.6 so 



that we find (0) = 0.021 {e'^/hnt) which is in good 
agreement with the computed gap. We find excellent 
agreement for other values oi d/t as well. This permits 
us to identify the higher-energy dispersion branch as the 
pseudospin-wave mode. Our animations''^ show that the 
other high-energy branch is an optical phonon mode in 
which the two bimerons in the unit cell have opposite 
motions. 
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FIG. 7: (Color online) Collective mode spectrum of a pseu- 
dospin meron crystal in a SLA configuration similar to that 
of Fig. 3(b) with v = 1.2, Asas/ {e^/^tj = and d/l = 0.8. 
The path in the Brillouin zone is as indicated in the inset of 
Fig. 4(b) and the symbols are as indicated in the legend of 
that figure. 

Figure 7 shows the collective mode spectrum of a 
pseudospin meron crystal. The two gapless modes are 
the phonon and the XY mode as discussed above and 
the other two branches are the optical phonon and 
pseudospin-wave modes. These two branches are degen- 
erate along TM at all d/^ for the SLA configuration. Eq. 
([55)1 fails in this case since (pi_3 (0)) = for a meron 
crystal (but (pi,3 (G ^ 0)) 7^ oj and we get (0) = 
for all interlayer separations. That this equation works 
for the bimeron crystal but fails for the meron crystal can 
be understood by the fact that the bimeron crystal has 
regions of coherent liquid between the bimerons whereas 
the meron crystal has no such regions. 



VII. CP^ SKYRMION CRYSTALS 

Now that we understand the limiting cases of pure spin 
and pure pseudospin crystals, we are ready to look at the 
collective excitations of the CP^ Skyrmion crystal with 
has entangled spin and pseudospin textures. We study 
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the excitation spectrum as a function of interlayer sep- 
aration, tunneling, and bias as was done for the ground 
state in Paper I^. We refer the reader to this paper for 
the definition of the various groundstates that we will use 
in the present section. We remark that all our data for 
the CP^ Skyrmion crystal in this section arc for filling 
factor ly ~ 0.8 where we get a good convergence of the 
order parameters for the CP^ phases. For > 1.0, a sim- 
ilar phase is difficult to find numerically (see Paper I^). 
For the spin-Skyrmion and pseudospin-Skyrmion crys- 
tals, our results above were for > 1.0. In these two 
cases, however, the phase diagram below = 1, can be 
related to that above v = 1 by particle-hole symmetry. 
In the CP^ Skyrmion case, there is no such particle-hole 
symmetry. 

The groundstate energy of Eq. ^ can be written in 
the spin-pseudospin langage as 

^ = -iP^ (0)) - ^ {S. (0)) - ^ (P. (0)) 



i- ^ Ti (G) |(p (G))|^ + i ^ J,,i (G) (G))|2 

G 

E ^(G)l(Sa(G))r (57) 

G a=R,L 

^E E ^(G)[l(^'.,a(G))|' + |(P,,„(G))|^ 
^X(G) [|(pm(G)>|V|(p2,3 (G))|^^ 



G a=^ 
1 



where we have defined the interactions 

Ti(G) = i/(G) + i/(G)-ix(G), 

J,,i(G) = H(G)-U{G)-\x{G)., 



and the fields (see Eq. 

1 



(G) 

P..- (G) 
Py.+ (G) = 
Py.- (G) = 



[(Pl.s) + (P3,l)] 



2 [(P2.4) + (P4,2)] 



2i 



[(Pl,3) - (P3,l)] . 



^ [(P2,4) - (P4.2)] 



In Eq. m,S, = Y.a=R..L S^,a and = Ea=- 
where i = x,y, z. The total density is defined by 



(58) 
(59) 

(60) 
(61) 
(62) 
(63) 

P 



{p (G)) = (pn (G)) + (p22 (G)) + (p33 (G)) -f (p44 (G)) . 

(64) 

We will also make used of the fields 

Rxx[c\) = ^ [pi,4(q) + P4a(q) + P2,3(q) + P3,2(q)I?5) 

^^^(q) = \ - P2,2(q) - P3,3(q) + P4,4(q)l66) 



A. Variation with interlayer separation 

At small value l^sAS I {e^ I ^^^) = 0.0002 << 
/S.Z / {g^ / 1^^) of the tunneling parameter, the ground 
state of the 2DEG at small interlayer separation d/t\s a. 
spin-polarized meron crystal with an SLA configuration^. 
The corresponding excitation spectrum is shown in Fig. 
8. For all the following dispersion plots, we represent 
the poles of Xnn by the full (black) circles, those of the 
pseudospin response functions xpp by the empty (blue) 
squares, those of the spin response functions xss by 
the right (red) triangles and by the left (green) trian- 
gles the poles of the two response functions Xj?a;x,-Rxi and 
XR^^,Rz^- The operator R^x (see Eq. (|55|)) involves a 
change in both the spin and pseudospin indices unlike the 
operators S, P or n. There are clearly two basic mode 
structures in Fig. 8. One, noted "(1)" is the pseudospin 
meron crystal response already studied in Fig. 7 and it 
comprises four branches. The second structure, noted 
"(2)", has two twice-degenerate branches for a total of 
eight branches for the whole CP^ crystal (considering the 
low-energy modes only since other higher-energy modes 
are also present in the GRPA). One branch in (2) is the 
Zeeman mode (gapped at Az/ {e^ / nt) = 0.004) while 
the other, which is also gapped, seems to appear pre- 
dominantly in Xfl^.-Rx. and XR..R.,- 
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FIG. 8: (Color online) Collective mode spectrum 
of a spin-polarized meron crystal with parameters: 
u = 0.8, Az/ (eV'i^) = 0.004, AsAs/ (eVK^) = 
0.0002, Ai,/ [e^/til) = 0, and d/l = 0.3. 

Figure 9 shows what happens to this collective mode 
spectrum when we increase the interlayer separation and 
keep both A5AS and fixed. In this case (see Fig. 3 
of Paper I^), the ground state becomes a CP^ crystal 
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at a critical layer separation dc that increases with Az- 
For Az/ {e^/nl) = 0.004, we find djt « 0.6. In Fig. 9, 
d/l — 1.0 and we see that the degeneracy of the modes 
in the structure (2) of Fig. 6 has been lifted and that two 
seemingly gapless modes have emerged from the lowest 
branch. In fact, the gap in the lowest branch in structure 
(2) vanishes continuously as the transition is approached. 
We believe that this should be observable experimentally. 
In the following analysis, it will become clear that of 
these two modes, one only one is truly gapless. We have 
identified the phonon, pseudospin-wave, pseudospin-yZ, 
and Zeeman modes in Fig. 9. We find that at d/l = 1.1, 
the modes of structure (2) become soft at point M in 
the Brillouin zone indicating an instability of the CP"^ 
crystal analogous to that of the uniform (v = 1) system 
at similar separations^. 
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and AsAs/ (e^/Kl) = 0.0044. From Fig. 5(c) of Paper I, 
the ground state of the 2DEG is a SS state in this case. 
The CP^ Skyrmion state is still stable even if is not the 
ground state and, as we see from Fig. 10, there is a very 
small gap with energy smaller than 0.001 (b^/k^) in one 
of the two modes of structure 2. As we will show below, 
the application of an electrical bias does not open new 
gaps in the system so that we conclude that with finite 
bias, Zeeman, and tunneling energies, our CP^ crystal 
has two Goldstone modes: the phonon mode and a phase 
mode whose origin we will discuss below. 
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FIG. 10: (Color online) Collective mode spectrum of a crys- 
tal of CP^ Skyrmions with strong tunneling. The param- 
eters are: u = 0.8, Az / {e'^ / k.£) = 0.008, Asas/ (eV«:^) = 
0.0044, Ai/ [e^/Kl) = 0, and d/i = 1.0. 



FIG. 9: (Color online) Collective mode spectrum of a 
crystal of CP'^ Skyrmions along path -M-X- in the Bril- 
louin zone. The parameters are: v — 0.8, Az/ (g^/k^) — 

0.004, AsAs/ (eV«:^) 
1.0. 



0.0002, At/ {e^/Ktj = 0, and d/i 



B. Variation with tunnel coupling 



The collective mode spectrum of a SS state aX d/l — 
1.0, Az/ (eV^^) = 0.008 and Asas/ (eV^f) = 0.0044 is 
shown in Fig. 11. For Az/ (c^/k^) = 0.008, the tran- 
sition to the SS state occurs at A^^^/ (e^//t£) w 0.004 



which is roughly A 



SAS 



0.5Az. The phonon and spin 



mode dispersions in Fig. 11 modes are typical of that 
of a Skyrme crystal (see Fig. 4(a)). There are only two 
gapless modes: the phonon and spin XY phase mode. 



To check if the two modes in Fig. 9 are really gapless, 
we plot in Fig. 10 the collective mode spectrum of the 
CP^ crystal at a higher value of tunneling. In Paper 
it was shown that when tunneling is increased, there is 
a transition from the CP^ crystal state to a symmetric 
Skyrmion (SS) state. The SS is basically a pseudospin 
polarized spin-Skyrmion state in which the 5, f and S, [ 
only are occupied. There is thus a spin texture on top 
of a state where all pseudospins point along the x axis. 
For Fig. 10, we choose d/i = 1.0, Az/ (eV^^) = 0.008 



C. Variation with electrical bias 

We now turn on the bias between the two wells so that 
all the charge is gradually transferred to the left well. 
If we start with a CP^ crystal at zero bias, our Hartree- 
Fock analysis (see Fig. 9 of Paper I^) shows that the spin 
and pseudospin textures are maintained until the charge 
has completely gone into one well and we get a spin- 
Skyrmion crystal. From the collective mode spectrum 
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FIG. 11: (Color online) Collective mode spectrum of 
a crystal of symmetric Skyrmions. The parameters 
are: v = 0.8, Az/ (eV«:^) = 0.008, Az/ (eV«:^) = 
0.006, At/ (eVw^) = 0, and d/^ = 1.0. 



FIG. 12: (Color online) Collective mode spectrum of a crys- 
tal of CP^ Skyrmions with a finite bias. The parame- 
ters are: v = 0.8, Az/ (eV«:^) = 0.01, Asas/ (eV«:^) = 
0.0002, At/ (eV«:^) = 0.20 and dji = 1.0. 



in Fig. 12, we see that the two Goldstone modes are 
preserved at finite bias if we consider that one of the low- 
energy mode has an extremely small gap (whose origin 
we explain in the next section). The qualitative features 
of the spectrum are unchanged from that of Fig. 9. The 
CP^ crystal is stable up to the largest bias for which 
there is a charge in the right well. When the charge is 
completely transferred to the left well, the spectrum is 
that of the Skyrme crystal shown in Fig. 4(a). 



where |0) is the vacuum state. We replace^ X by r in 
this last equation so that our order parameters can be 
written in real space as 



D. Physical origin of the Goldstone modes 



(P<Ti,^2 W) = <, (r) fl'T^ (r) 



(69) 



Following Gosh and Rajaraman^, we pametrize a gen- 
eral CP^ spinor by six angles that can be roughly inter- 
preted as the polar and azimuthal angles of the spins in 
the right and left wells: du^^h and yj^j, (^l, and the polar 
and azimuthal angles of the total pseudospin, a and /?. 
Writing 



We have 



cos § cos 
COS f sin ^-fe^VR 
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(67) 



a possible spin-pseudospin texture state can be written 
as 



X 



'<y,X 



|0), 



(68) 
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and 
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(Pi,3 (r 
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(P2,4 (r 
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(71) 
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and all angles are functions of r. 

The symmetry-breaking fields in our Hamiltonian in 
Eq. ^ couple to Px, Pz, and 5*2 where 



{S. (r)) = - 
{P. (r)) 



sm — cestui + cos —cos On 



(72) 



■ sm a 



- sma 
2 



cos COS COS p 

2 2 



sm — sm — COS (^ + - (^i?J 



(^'^ W) = -cosa. 

These expressions may be substituted into our Hartree- 
Fock Hamiltonian of Eq. ([2]). The origin of the second 
Goldstone mode (the phase mode) is then clear: the en- 
ergy of the CP^ crystal is unchanged if we change the 
azimuthal angles (p/? and by the same amount (i.e., 
keep ipLt — <pl fixed) for all spins, whereas, in the CP^ 
structure, the symmetry in Lp^i -I- 1^9^ is broken. The sec- 
ond Goldstone mode of the CP^ crystal is thus an XY 
phase mode corresponding to a global rotation of the 
spins at q = i.e. an XY phase mode in Lp^t + ipL- 

Were it not for the term cos {(3 + lpl — fn) in the ex- 
pression of Px in Eq. (|72|) , our Hamiltonian would also be 
invariant with respect to ipfi — ipL and to /3: we could then 
have two additional Goldstone modes. The tunnel cou- 
pling is of the form As as J dv (P^ (r)). In our numerical 
solutions we find that the CP^ crystal phase is stable at 
very small values of Ag^s and, furthermore, at these val- 
ues (Px (r)) is also small. For example, in the case of Fig. 
12, we have / dr {Px (r)) = 0.005 and Asas/ (eV^^) = 
0.0002 so that a change in (pn — (pL or f3 would lead to a 
change in the energy (a bare gap) Asas J dr {Px (r)} by 
a quantity smaller than 1 x 10~® [e^ / k,£^ which is smaller 
than our numerical accuracy. Even if we allow for this 
change in energy to be renormalized by vertex correc- 
tions (as is the case for the pseudospin-wave mode in the 
UCS in Sec. IV where the gap AJ^g ~ lOAgyfis), the ex- 
pected gap would still be small. This is the reason why 
there seems to be three gapless modes in Fig. 9. In Fig. 



10, we have choosen Asas I (e^//t£) = 0.0044 and obtain 
/ dr {Px (r)) = 0.163 in this case. We could thus expect a 
bare gap of order 0.0007 {e^ / nt) which is consistent with 
the gap we find numerically. We thus associate the small 
gap mode with either ipn — ^pl or /3 as these angles are 
present in the expression of both S_l and Pj^. Presum- 
ably one of these angles (or a combination of them) leads 
to the mode with the small gap 

One last point is to check that a change ^pr,^l or /3 
does not change the Skyrmion density. We recall that the 
topological charge density of CP^ skyrmion is defined^ 
by 



q{v) 



27r 



(73) 



where is the covariant derivative of the U{1) gauge 
transformation with the gauge defined by (r) = 
a*d^afj i.e. — d^-\- iA^. With the parametriza- 
tion given by Eq. ([57]) . we have 
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TT ^ — ' 



(74) 



which leads to 
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The topological density is related to the charge density^ 
so that a global change in the azimuthal angles has no 
effect on the charge density itself. 



VIII. CONCLUSION 

In a previous work^, it was shown that a crystal with 
entangled spin and pseudospin textures, i.e. a CP"^ 
Skyrmion crystal, can be the Hartree-Fock ground state 
of the 2DEG in a DQWS when interlayer coherence is 
established at small interlayer separation. In the present 
work, we have presented an exhaustive study of the col- 
lective excitations of such a crystal by working in the 
generalized Random-Phase approximation. In the limit- 
ing cases of a pure spin-Skyrmion crystal i.e. at strong 
tunneling or at strong bias, we found that the spin- 
Skyrmion crystal has a gapless phonon mode and a sep- 
arate Goldstone mode that arises from a broken U (1) 
symmetry. At zero Zeeman coupling, we demonstrated 
that the constituent Skyrmions broke up, and the result- 
ing state is a meron crystal with four gapless modes. In 
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contrast, a pure pseudospin Skyrme crystal at finite tun- 
neling has only the phonon mode as a Goldstone mode. 
For A5yi5 0, however, the pseudospin Skyrme crys- 
tal evolves into a meron crystal and it then supports an 
extra gapless {U{1)) mode in addition to the phonon. 

For a CP^ Skyrmion crystal, we found a U{1) gapless 
mode (in addition to the phonon mode) in the presence 
of non-vanishing symmetry-breaking fields As^s (tun- 
nehng), Az (Zeeman coupling), and (electrical bias). 
This phase mode is explained by the fact that the en- 
ergy of the CP^ crystal is unchanged if we change the 
azimuthal angles tpu and fL oi all the spins in the right 
or left layers by the same amount. In addition to this 
gapless mode, we found a second mode with a very small 
gap. We believe that this mode can be associated with 
the variation of the relative azimuthal angles (pu — (pL ot 
with the variation of the azimuthal angle P of the total 
pseudospin vector or to a combination of these angles. 

Having established the low-energy excitations of the 



CP^ Skyrmion crystal and obtained the corresponding 
response functions, it is now possible to compute the 
NMR relaxation time Ti and see how it compares with 
the measured experimental values^ in order to see if the 
formation of a CP^ crystal plays a role in these experi- 
ments. This, however, is beyond the scope of the present 
work. 
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